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Abstract 

Abstract Low-T Brillouin spectra of the incipient ferroelectric KTaOs 

I exhibit a broad central peak (CP), and additional Brillouin doublets (BD), 

^ , that can both be related to phonon-density fluctuations. On the basis of 

Q ■ new high-resolution neutron data obtained of low-lying phonon branches, 

O \ we analysed the phonon-kinetics mechanisms that are possibly the origin 

of these unusual features. Firstly, transverse acoustic (TA) phonons whose 

normal damping is faster than the BD frequency can produce hydrody- 

namic second sound. Secondly, two-phonon difference scattering from 

low damping thermal transverse phonons contribute to the spectra with 

either a sharp or a broader doublet, depending on the phonon group ve- 

Qv \ locity and anisotropy of dispersion surfaces. The position of the observed 

C^ • sharp doublet is consistent with both mechanisms, but a comparison of the 

^^ I computed and experimental anisotropics favours the second process. 

o 

"^ I Keywords quantum paraelectrics, inelastic neutron and Brillouin scat- 

tering, second sound, two-phonon difference scattering, phonon-density- 
fluctuation processes. 
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Both KTaOs and SrTiOs perovskite-type crystals exhibit ferroelectric transition- 
like features with decreasing temperature. Their dielectric constant, e, increases 



;_j ■ as the energy Usm of the associated Brillouin zone (BZ) centre transverse optic 



(TO) mode decreases following the classical Curie-Weiss law. However, zero- 
point quantum fluctuations prevent critical mode condensation and the crystal 
remains paraelectric down to the lowest temperatures (T). Such a material is 
then called a quantum paraelectric (QPE) [[I], ^ ^]. The proximity of the soft TO 
and acoustic modes near BZ centre might significantly modify low-T phonon 
kinetics. The case of SrTiOs is more complicated than KTaOs, since SrTiOa 
undergoes a well-known structural transition at Tq ~ 105 K [Q,!], ^. Below this 
temperature, the unit cell doubles and the crystal usually forms domains. In ad- 
dition to this experimental difficulty, the description becomes more involved due 
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2 THE BROAD BRILLOUIN DOUBLET AND CP OF KTaOs 

to the presence of three additional low frequency branches. For this reason, we 
first investigate KTaOs, which remains in a cubic PtuSui space group symmetry 
down to the lowest temperatures. 

Many inelastic neution-scattering [^ pi j|], Raman jj^], infrared [11, [l^ 



Rayleigh-Brillouin [13, |l^, 15], dielectric []T6], thermal [ jT?! ] or mechanical [ |18| ] 
studies highlighted numerous anomalies in pure KTaOa crystals at low temper- 



ature. Particularly, Hehlen et al. [14, 15] observed an extra Brillouin doublet 



(BD) with a sound-like dispersion relation in KTaOs and SrTiOa. In KTaOa, it 
is seen for T < 25 K along <001> and <110> directions. This new excitation 
seems to appear on top of the quasi-elastic broad central peak (CP) first reported 



by Lyons and Fleury [ ]13| ] 



The aim of this paper is to analyse the BD and CP anomalies in pure KTaOs 
in terms of classical phonon kinetics, /. e. in terms of phonon-density-fluctuation- 
induced light-scattering. In Section ^, we present some new Brillouin measure- 
ments studying the anisotropy of the CP and doublets along the three principal 
symmetry axes. High-resolution data for low energy dispersion branches are ob- 
tained using inelastic neutron-scattering, as described in Section^ Experimental 
results are then used to describe low energy phonon dispersion surfaces around 
the BZ centre by mean of a simple phenomenological parameterisation. Using 
such a model, we shall show (Section ^ that the anomalous Brillouin contri- 
butions may originate from either two-phonon difference scattering (TPDS) or 
second sound (2S) processes, but the observed anisotropy seems to favour the 
former. 



2 BRILLOUIN SCATTERING EXPERIMENTS 

The KTaOs Ultra High Purity crystals were grown by one of us (LAB) at ORNL 
by 'modified spontaneous nucleation' within a slowly cooled flux. The light- 
scattering spectra were excited using a single-mode argon-ion laser operating 
at 5 145 A and observed in backscattering configuration along the three crystal 
principal directions. The momentum exchange was \q\ = 5.6 10^'^ A^^. The 
temperature T was varied between 5 and 300 K. The spectra were analysed by 
mean of a six-pass tandem interferometer that provides an excellent contrast and 



suppresses the overlap of orders [19] 



Figure |l| shows typical spectra for various temperatures and q// <110>. For 
T > 100 K, an intense CP appears, strongly coupled with the LA phonon (curve 
A). As the temperature is decreased, the CP contribution becomes weaker, and an 
unexpected new broad Brillouin doublet (BBD) grows upon its shoulders (curve 
B). For T < 25 K, the sound-like BD contribution, first observed by Hehlen et 



al. [ ]14| ] below 22 K, can be measured (curve C), with a corresponding speed 



vbd ~ 1 100 m/s. It exhibits a width of typically Tbd ~ 6 GHz. Along the 
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Figure 1 : Brillouin backscattering spectra of KTaOa showing the acoustic LA 
phonon, central peak (CP), Brillouin Doublet (BD) and Broad Brillouin Dou- 
blet (BBD), measured at (A) 195 K, (B) 31.5 K, and (C) 13 K with \q\ II <110> 
principal axis. Empty zones correspond to central laser peak and first order in- 
strumental ghosts [|19|]. Spectra are shifted for clarity. Along <001> and <111>, 
one can observe shapes A-C and A-B respectively. 



<001> direction, the BBD excitation is not seen, while along <111>, no BD 
appears at temperatures as low as 5 K, and the CP is broader than along other 



directions below 100 K []20|]. The CP is clearly coupled to acoustic phonons 
along the <110> and <111> directions. 



3 LOW ENERGY PHONON DISPERSION CURVES 
AND THEIR PARAMETERISATION 

In order to analyse the origin and the anisotropy of excitations CP, BD and 
BBD in terms of phonon based processes, we parameterised the low energy 
phonon dispersion surfaces, which are still significantly populated under 100 
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Figure 2: Dispersion curves of KTaOs low energy phonons between 10 and 80 
K along high symmetry axes. The mean energy uncertainty is 0.1 meV. The 
LA mode is for <0 2+^>. Lines correspond to calculations using the Vaks 
parameterisation. 



K in KTaOa. A more ambitious modelisation has already been performed along 
high symmetry axis [^ using the polarisability model [^, ^] but it requires both 
substantial computations and the determination of 15 parameters in order to be 
applied just along the three principal directions. 

As the precision of the previous studies were not sufficient for our purpose, 
we performed new high-resolution inelastic neutron-scattering measurements on 
IN 14 at the Institut Laue Langevin (ILL), Grenoble, France [2^, 22]. By taking 
into account a 4-dimensional modelisation of the instrumental resolution ellip- 
soid and of the anisotropic phonon dispersion surfaces, it was possible to mea- 
sure phonon branches (Figure ^ with a typical energy resolution of 0. 1 meV. 
This enabled a direct comparison of the soft mode energy (2.6 meV at 10 K) and 



damping (0.07 meV at 10 K) with Hyper Raman results near BZ centre [10]. 
Some off-principal axis experiments were also performed. Acoustic phonon 
dampings could not be measured, setting a maximum value of about 0.05 meV 
from instrument resolution. 

These Brillouin and neutron-scattering measurements, enabled us to apply 
the Vaks parameterisation |23, 0, 25] on the five lowest phonon modes for tem- 
peratures below 100 K, in the whole central part of BZ {\q\ < 0.2 rlu) where 
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phonon dispersion surfaces are significantly populated ||20|, g2|]. Such a model 
requires the determination of only 8 parameters, 4 of which can be measured 
directly by mean of light-scattering. The resulting curves along high symme- 
try axes are shown in Figure ^. As presented in Figure ^-a, transverse phonon 
dispersion surfaces are quite anisotropic, and show deep 'valleys' (low energy 
and group velocities), induced by TA-TO coupling [^, and 'uphill' (high energy 
and group velocities) regions depending on mode propagation and polarisation 
directions. 



4 PHONON-DENSITY-FLUCTUATION PROCESSES 



Hehlen et al. [15] considered several possible explanations for the origin of 
the new BD. The most likely one at the time seemed to be the onset of second 
sound (2S), a phenomenon independently predicted in QPE by Gurevich and 
Tagantsev [|2^]. The physical reason for this effect is that the presence of the low 
frequency TO branches enhances normal (momentum conservation) relaxation 
processes. This could open the frequency 'window' condition F/j < ij02S < 



Ftv [15, 26 1 for the appearance of 2S near BZ centre, where V^ and Vr are 
the normal and resistive dampings respectively. In KTaOa, this excitation is 
expected to be nearly isotropic in velocity and to couple to acoustic modes. 

Following a thermodynamical approach, it is possible to evaluate the 2S ve- 
locity near BZ centre using t;|g = -^^, where S is the entropy density, C is 



the heat capacity, and P is a density tensor, as detailed in [26|. Phonons con- 
tributing to 2S should have normal dampings higher than the doublet frequency 
ojbd, in order that the 'window' condition is satisfied. It is possible to evaluate 



three-phonon normal dampings [27| through electrostrictive processes [18, 24 1 
using the phonon dispersion surfaces parameterisation. Numerical computations 
of this sort indicate that TA dampings are minimal in <001> and <1 10> 'valley' 



directions, and higher elsewhere, particularly around 'uphill' regions [[20|]. 

Computations of 2S velocity using the Vaks model and the thermalized part 
of the phonon spectrum give V2s = 1 100 ± 200 m/s, in excellent agreement with 
Brillouin scattering measurements of BD speed. However, it seems difficult to 
explain the directional dependence of BD (see Section ^ and Figure [T]) in terms 
of2S. 

Second sound is the collective mode of phonon density fluctuation processes 



[28]. When the window condition is not fulfilled {lo2S > ^n), non-thermalized 



phonon can still interact via two-phonon difference scattering (TPDS) processes. 



Following Wehner and Klein [|28|], the resulting Brillouin second order scattering 
intensity is given by [EO] : 
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Figure 3: (a) Low energy phonon spectrum at T = 10 K from Vaks parameter- 
isation for a rotation |^| = 0.05 rlu. Directions <110>, <111>, and <001> are 
for angles 0° or 180°, 35.5° or 144.8°, and 90° respectively. Off principal-axis 
neutron-scattering data are also indicated, (b) Computed two-phonon difference 
scattering spectrum from different TA model dispersion sheet shapes at T = 23 
K. Phonon damping is 1 GHz . 
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where n is and 4 for TO and TA contributions respectively. This expression 
has been evaluated for model phonon dispersion surfaces in the case of valley, 
isotrope and uphill shapes (see Figure ^b). It is then possible, in the case of 
low phonon dampings, to obtain a sharp peak from valley contribution, or a 
broader quasi-central spectrum from uphill contributions. The peak or thresh- 
old frequency u>tpds is fixed by the mean group velocity of contributing non- 
thermalized phonons, and is smaller in valley than uphill. The TPDS spectrum 
broadens with increasing dampings, and then appears as a broad doublet or cen- 
tral peak, depending on group velocity distribution. Thus, the TPDS spectrum is 
essentially anisotropic. 

In this scheme, the BD low temperature contribution originates from low 
damping and group velocity 'valley' TA phonons TPDS i.e. for <001> and 
<110> propagations (see Figure ^a). The BBD comes from higher damping 
'uphill' TA phonons TPDS i.e. for <110> and <111> propagations, at slightly 
higher temperatures where those regions are significantly populated. The BD and 
BBD contributions disappear when the soft mode energy increases, and the TA- 
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TO coupling induced valleys open up. For the TO modes, with high dampings, a 

very broad quasi-central peak is expected, corresponding to the CP contribution. 

Wehner and Klein [28] predicted that in various crystals, including SrTiOs, 



the integrated intensities of phonon density fluctuation processes (i.e. 2S or 
TPDS) could be 10 to 100 times the value given by the Landau-Placzek ratio. It 
is then quite plausible that these intensities become comparable, or even higher 
than phonon light scattering intensities in KTaOa, which would thus bring them 
in agreement with observed intensities. 

5 CONCLUSION 

We have shown that the phonon spectrum and dampings in KTaOa are consistent 
with two phonon-density-related scenarios for the appearance of quasi-elastic 
central peak and doublets. One scenario is related to temperature waves (second 
sound) involving a part of the spectrum that can be thermalized during the dou- 



blet lifetime [15, 26]. The other scenario is related to the two-phonon difference 



scattering of non-thermalized transverse phonons from phonon sheet 'valleys' or 



'uphill' regions [2C]. The position of the sharp doublet BD is consistent with 



both scenarios whereas the anisotropy of the scattering favours the second one. 
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